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We calculate the contribution of cosmic strings arising from a phase transition in the early universe, 
or cosmic superstrings arising from brane inflation, to the cosmic 21 cm power spectrum at redshifts 
z > 30. Future experiments can exploit this effect to constrain the cosmic string tension G/i and 
probe virtually the entire brane inflation model space allowed by current observations. Although 
current experiments with a collecting area of ~ 1 km'^ will not provide any useful constraints, future 
experiments with a collecting area of 10" - 10^ W covering the cleanest 10% of the sky can in 
principle constrain cosmic strings with tension G/i ^ 10"^" — lO"^'^ (superstring/phase transition 
mass scale > 10^^ GeV). 



Introduction. — There has been a revival of interest in 
cosmic strings, the line like topological defects of cosmic 
length, after it was found that they can arise in the su- 
perstring theories in braneworld inflation scenarios [ij. 
They are also found to be inevitable in a wide class of 
supersymmetric Grand Unified Theories (GUTs) 0. His- 
torically cosmic strings were found to form in GUTs dur- 
ing phase transitions in the early Universe along with 
the other topological defects like monopoles and domain 
walls 0]. Unlike monopoles and domain walls, which 
very quickly dominate the energy density of the Universe 
if formed after inflation, strings ap;^oach a scaling solu- 
tion and can remain sub-dominant 4] . In brane inflation 
only cosmic strings and no monopoles or domain walls 
are produced @. 

Cosmic strings were proposed as a mechanism for gen- 
erating the primordial fluctuations which later grew to 
form the large scale structures we see today Q- The 
Cosmic Microwave Background (CMB) anisotropics and 
the matter power spectrum arising from the fluctuations 
seeded by cosmic strings are very different from those 
generated from inflation fT*]. Inflation just prescribes the 
initial fluctuations at the end of inflation generated once 
and for all, which then just evolve. Cosmic strings gener- 
ate fluctuations throughout the history of the Universe. 
One effect of this is that the fluctuations generated at dif- 
ferent times add up out of phase and wash out the acous- 
tic oscillations in the CMB. The discovery of the acoustic 
peaks by CMB experiments d, Q was a major success 
for inflation and ruled out cosmic strings as the dom- 
inant mechanism for seeding the cosmic perturbations. 
A sub-dominant contribution from cosmic strings to the 
cosmic perturbations is still not ruled out with the cur- 
rent constraint being G/i < 10~^ for classical strings [lo| . 
Cosmic strings if discovered either through their gravita- 
tional lensing effects [ll| , through the gravitational waves 
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produced at string cusps or decaying loops 12] or through 
their effect on the CMB and the matter power spectrums 
will provide insight into the fundamental physics at high 
energies which is beyond the reach of currently planned 
terrestrial experiments. 

21 cm radiation from z > 30 is an excellent probe of 
the state of the Universe at that time. This radiation 
can probe much smaller scales than the CMB, in the 
redshift range 30 < z < 200, and provides a three di- 
mensional view of the Universe before reionization p^ . 
It has been shown to be an excellent probe of the fun- 
damental physics like variation of constants (fractional 
variation in the flne structure constant of < 10^"" with 
10** km^ collecting area), non-Gaussianity from inflation 
(non-Gaussianity parameter /„; 0.01 — 1), dark mat- 
ter and inflationary fundamental physics [ij, [3, [IH . In 
this Letter we show that it can, in principle, put unprece- 
dented tight constraints on the cosmic string contribution 
to the perturbations in the matter or equivalently on the 
string tension G/i, and possibly other parameters, which 
translates into constraints on the GUTs and the super- 
string theory. G is the gravitational constant and /i is the 
string mass per unit length so that G/i is dimensionless. 

Cosmic strings. — Cosmic strings arise in GUTs and su- 
perstring theories whenever there is a phase transition in 
the Universe if the vacuum manifold contains unshrink- 
able loops, e.g. U(l). The superstring theory can pro- 
duce a variety of cosmic strings, which can be fundamen- 
tal (F-)strings or D-strings produced during annihilation 
of D-branes. The string tension for these st ring s in brane 
inflation models is lO'^^ < Gfi < 10"^ Sill [11. Just 
like the classical cosmic strings from GUTs, they are wig- 
gly, can intercommute and form loops which can decay 
into gravitational radiation or elementary particles. The 
main difference from the classical cosmic strings of GUTs 
is that their intercommuting probability can be less than 
unity. Also different kinds of strings can form bound 
states 119^. The superstring theory string networks have 
scaling solutions 0| just like the classical strings ^ i.e. 
the total length of the strings inside a horizon volume 
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is proportional to the horizon size at any time and the 
string energy density is a constant fraction of the dom- 
inant energy density component of the Universe in the 
radiation dominated as well as matter dominated eras, 
preventing them from dominating. This also makes possi- 
ble to construct simpler models of string networks which 
can than be used to study their impact on cosmology. 
Also since the perturbations produced by cosmic strings 
are independent of the inflationary initial conditions, the 
two kind of perturbations can be evolved independently 
and the resulting power spectra for CMB or matter added 
together to get the total power spectrum. 

The long wiggly cosmic strings have a structure that 
resembles a random walk on scales lar ger than the hori- 
zon but are straight on small scales [2l|. We use the 
CMBACT code developed by Pogosian and Vachaspati 
[13] which is based on CMBFAST [H]. Wiggly strings 
are modeled in the code as independent pieces of small 
strings whose length is taken to be of the order of the 
correlation length of the pieces of strings derived from 
numerical simulations. The intercommuting probability 
(P) is taken to be unity. The wiggliness on small scales 
results in the effective string tension and mass per unit 
len gth of the string observed by a distant observer to dif- 
fer [2J| with the equation of state UT = ^"^ , where IJ and 
T are the effective mass per unit length and tension of 
the wiggly string. Thus at scales greater than the scale 
of wiggles, the matter around experiences a Newtonian 
gravitational potential in addition to the deflection due 
to the conical space around the string [IHl- The wiggli- 
ness can be controlled in CMBACT by a wiggliness factor 
a defined by the equations IJ = a/i and T = ^/a (2^ . 
21 cm radiation- After recombination most of the bary- 
onic matter is in the form of neutral atomic hydrogen and 
helium with a very small residual ionization. This small 
fraction of electrons couples the CMB to the matter up 
to z ~ 500 through Compton scattering and maintains 
the gas at the same temperature as the CMB. Around 
z ~ 500 the Compton scattering timescale becomes larger 
than the Hubble time and the gas decouples from the 
CMB and cools adiabatically thereafter. The ground 
state of hydrogen atom has a hyperfine splitting with 
an energy difference of — 0.068 K. This corresponds 
to the V21 1420 MHz rest frequency or A21 ~ 21 
cm rest wavelength. The hydrogen gas will absorb or 
emit photons of this energy depending on the popula- 
tion levels in the two states which is best parameter- 
ized by the spin temperature Tg defined by the equation 
TLt/ris = gtl gs^~'^*^^' ■ Here nt and Us are the number 
densities of hydrogen atoms in the excited triplet state 
and the ground singlet state respectively, gt = i and 
— \ are the corresponding statistical weights. The 
population levels during z > 30 can change through col- 
lisions or through emission and absorption of CMB pho- 
tons. Initially the collisions dominate over the radiative 
process and follows the gas temperature Tg. At late 
times the density of gas becomes too low for collisions 
to be effective and approaches the CMB temperature, 



T~f = 2.725(1 -I- z). Thus in the redshift range of about 
500 > z > 30, Ts < T-y and we have a net absorption of 
the CMB [H,!!!]. At z 30 the first stars are born and 
the evolution of the Universe enters non-linear regime. 
We will focus on the redshift range 200 > z > 30 in this 
Letter, which corresponds to the observed frequency of 
- 7.1MHz <v < 47.3MHz. 

The observed intensity ly can be expressed in terms 
of the brightness temperature using the Rayleigh-Jeans 
formula, = l^c? j^ksv^ , with T;, given by [13] 

^ _ (Ts - T^) T 3c^?tAiong 

(1 + ^) ' ^~ lQkBvUH+{l + z)f^)T; 

where nu is the number density of atomic Hydrogen, Aiq 
is the Einstein A coefficient for spontaneous emission, c 
is the speed of light in vacuum, ks is Boltzmann's con- 
stant and fi = h/2TT with h being Planck's constant. H 
is the Hubble parameter, r is the comoving distance to 
redshift z and v is the peculiar velocity along the line of 
sight. We ignore the contribution of the vector modes to 
peculiar velocities since it is sub-dominant (< 1%) com- 
pared to the scalar mode contribution at scales of interest 
{I > 1000). All quantities except the physical and atomic 
constants are functions of z. There will be spatial fluctu- 
ations in T{, and caused by the fluctuations in uh and 
Tg, which in the redshift range of interest are related to 
the linearly evolved primordial perturbations in standard 
inflationary cosmology. We will see below that the cosmic 
strings, if they exist, can add a significant contribution 
to these fluctuations. 

Expanding the fluctuations in the brightness temper- 
ature 6T1, — Th — Th, where Tf, is the mean brightness 
temperature, in spherical harmonics, we get the angu- 
lar power spectrum C;(z) = (a/ma/*„), where aim are the 
coefficients in the spherical harmonic expansion of STh. 
Following [1^ we can write, 

P^{k,z)Si{k,z), 

(zTTj 

where P{k, z) is the baryon (Fourier) power spectrum, 
index i — ad or cs for adiabatic perturbations from in- 
flation or perturbations from cosmic strings respectively. 
We have incorporated the 21 cm physics into Si{k, z) [2q . 

Results- We calculate P'"^{k,z) using CMBFAST [H 
and P^'ik, z) using CMBACT The cosmological pa- 
rameters are from WMAP3 assuming KG DM cosmology 
[§] . For the cosmic string model in CMBACT we use ini- 
tial rms velocity of 0.65, wiggliness factor in the radiation 
era of 1.9 and initial correlation length of 0.13 times the 
initial conformal time, motivated by numerical simula- 
tions [m, dl] . The intercommuting probability of strings 
is taken to be unity, which means classical strings. The 
effect of smaller intercommuting probability will trans- 
late into a denser network which will make G/i smaller 
for the same amplitude of string generated perturbations. 
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FIG. 1: Angular power spectra from inflationary adiabatic 
initial conditions (COBE normalized) and cosmic strings. 
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FIG. 2: Constraints from current and future experiments on 
Gjj for a sky fraction of 10%, bandwidth of 0.4 MHz and inte- 
gration time of 3 years. Also marked is the constraint on Gfi 
achievable by a full sky cosmic variance limited experiment 
for the parameters assumed in the present calculation. 

Our predicted constraints on Gfi are therefore conserva- 
tive. 

The angular power spectra for two redshifts are pfotted 
in Fig. [l]for G/i = 10~^ afong with the inflationary spec- 
tra. One important feature of the cosmic string power 
spectrum is that it turns over at smaller scales compared 
to the inflationary adiabatic power spectrum. This is due 
to the fact that the strings continue to generate pertur- 
bations actively at all times. To estimate the constraints 
possible on Gfi from these observations we calculate the 
Fisher matrix Fg for the parameter 9 = (G^/IO^^)^, as- 
suming that all the other parameters are known. 

z 1=1 ^ ^ 

where C; — Cf^ + Gf*, fsky is the sky fraction observed 
and the sum is over all redshift slices and / up to a max- 
imum I = i. We use the fact that the bandwidth of any 
experiment is finite only to estimate the noise Cj^ and the 
number of redshift slices available and ignore its damping 
effect on G/ [29]. Bandwidth of 0.4 MHz ensures that the 
slices are seperated by ^ r/l for I > 1000 and are thus 
uncorrelated. 



The noise in the frequency range 1 < v < A7 MHz is 
dominated by the sky temperature which follows a power 
law, (G/^)i/2 oc Tsky oc v-'^tJM- For LOFAR [sH we 
use the noise estimates of [32] at z=10, scaled to low 
frequencies using the above mentioned power law with a 
bandwidth Ai^ — 0.4MHz and inte grat ion time of 3 years. 
These results also apply to LWA [33'] which has similar 
specifications as LOFAR. Fig. [2] shows the constraints 
on assuming an error of l/2(i^g)^/^ (the factor of 2 
comes in when we convert the error on 9 to the error 
on 0^/^) for LOFAR, two futuristic experiments and the 
cosmic variance limit assuming that the foregrounds can 
be removed at required precision [s^, HI] . Most of the in- 
formation on G/i is at / > 10^ as is clear from an inspec- 
tion of the power spectra also. For the cosmic variance 
limit we use the fact that there are more independent 
modes at high 1, « v^rH / V2ilc [l^ to calculate the 
number of redshift slices. The curve labeled (100 km)^ 
corresponds to a futuristic telescope of size 100 km and 
collecting area of 10* km^ that will reach out to I ^ lO'^ 
and GfjL 10~^°. A (1000 km)^ telescope will be needed 
to constrain G^i ^ 10~^^. Such a telescope might be pos- 
sible not only on Earth but also in space [35| ] or on the 
far side of Moon [3^. To reach G^ ^ 10" will require a 
collecting area of 10^^ km^. Reaching the cosmic variance 
limit of ^ 10~^^ may not be possible because the scat- 
tering of radio waves by the ionized interstellar medium 
will limit the smallest angular scales (oc v^"^) that can be 
observed [13] • In particular the I > 10® modes will not 
be available at all redshifts. 

It is clear from Fig. [2] that the information content of 
the 21 cm signal is huge and can in principle constrain 
G/i ^ 10"^®, if we have a cosmic variance limited mea- 
surement on the full sky. This corresponds to a phase 
transition energy scale of 10^^ GeV for GUT theories. If 
we take /x ~ 2M^ for D-brane strings, this means bounds 
on the superstring mass scale, Mg, down to ~ 10^^ GeV 
[^. In reality only a fraction of the sky would be avail- 
able due to our being confined to the galaxy and small 
scale modes may not be available because of the inter- 
stellar scattering of radio waves, but even then the 21 
cm signal has impressive constraining power over Gfi. 
The power at small scales due to cosmic strings is in fact 
underestimated in this linear calculation. Cosmic strings 
generate wakes behind them, which have a density con- 
trast of unity, that is not taken into account here and 
which would enhance the power due to cosmic strings at 
small scales through non-linear gravitational evolution. 
This is just the information contained in the power spec- 
trum. Higher order correlations will provide additional 
discriminating power to check for the signatures of per- 
turbations seeded by a network of cosmic strings. The 
non-Gaussianities due to cosmic strings would be larger 
due to the highly non-linear nature of the perturbations 
at small scales and different from those produced during 
inflation [l^ . This impressive constraint is due the fact 
that there are large number of modes available at high 1 
so that the statistical errors become very small. 
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High redshift 21 cm observations thus provide a rare 
observational window into the superstring theory and 
supersymmetric grand unified theories through cosmic 
strings. 
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